Disorder, Order, and Domain Wall Roughening in the 2d Random Field Ising Model 
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Ground states and domain walls are investigated with exact combinatorial optimization in two- 
dimensional random field Ising magnets. The ground states break into domains above a length 
scale that depends exponentially on the random field strength squared. For weak disorder, this 
paramagnetic structure has remnant long-range order of the percolation type. The domain walls are 
super-rough in ordered systems with a roughness exponent C, close to 6/5. The interfaces exhibit rare 
fluctuations and multiscaling reminiscent of some models of kinetic roughening and hydrodynamic 
turbulence. 
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The random field Ising model attracts interest since it 
presents an example of competing mechanisms for order 
and disorder. The local spin couplings favor ferromag- 
netic ordering whereas variations in the random fields 
favor disorder. This competition affects thermodynamic 
properties. Disordered systems are very often governed 
by the zero-temperature behavior, or the structure and 
energy of the ground state (GS). In the random field 
Ising model (RFIM) statistical mechanics of interfaces or 
domain walls becomes the key question. This is impor- 
tant for finite-temperature dynamics, for coarsening, ag- 
ing, and fluctuations. The interest is in the scaling and 
universality in a problem dominated by a complicated, 
multivalley energy landscape [|l]||] . 

The effect of dimensionality on order or disorder was 
solved by Aizenman and Wehr. They proved rigorously 
that the two-dimensional (2D) RFIM has in the ther- 
modynamic limit no long-range ferromagnetic character 
||. In 3D order was shown to exist at finite tempera- 
tures and weak fields B , and therefore the lower critical 
dimension of the RFIM is 2. Here we study the tran- 
sition from order to disorder in the 2D RFIM: ground 
states and scaling properties of single domain walls with 
varying system size. The expectation is that the ground 
state becomes unstable to domain formation at a breakup 
length scale |i| because of domain wall entropy even at 
zero temperature. In contrast, a simple energy or Imry- 
Ma domain argument indicates for weak fields in two 
dimensions long-range order (if hL d / 2 < JL d ~ l , where h 
is the RF strength, J the strength of ferromagnetic cou- 
plings, L length scale, and d the dimension) This 
fails so that domains, albeit large ones, do exist for ar- 
bitrarily weak fields. Existing finite temperature Monte 
Carlo and exact ground state results || do not extend 
into this regime. 

The scaling properties of domain walls is studied here 
with the domain wall renormalization group (DWRG). 
One considers DW's imposed with boundary conditions, 
and compares to systems without forced DW's and with 
the same disorder to find the DW energy. The do- 
main walls are predicted to be self-affine by functional 
renormalization group calculations and an Imry-Ma ar- 
gument, with the roughness exponent C,^ F = (5 — d)/3 
P,^0|. C,f F shows, e.g., how the interface width scales, 
w ~ Z> (w 2 = (z 2 — z 2 ), where z is the local interface 



height). In this picture C vanishes at the upper criti- 
cal dimension and at d = 2 (2 = 1 The domain 
wall energy would concurrently be expected to be linear, 
E(L) = EqL + E\L 6 . The energy fluctuation exponent 6 
should obey the exponent relation = 2( + d — 3, sim- 
ilar to the random bond Ising model and directed poly- 
mers pT|,p|. The 1+1 -dimensional RF domain wall prob- 
lem maps in the continuum limit directly to the Kardar- 
Parisi-Zhang or Burger's equation, the paradigm of inter- 
face models in disordered media |i~2| , so that again (2 = 1 
and 2 = 1 @. 

The picture of self-affine DW's has been claimed to be 
confirmed by both early transfer-matrix calculations |l4| ] 
and studies using combinatorial optimization |l^,[l5| . In 
this Rapid Communication this is shown to be false. The 
domain walls exhibit rich scaling reminiscent of turbulent 
behavior as in certain kinetic growth models fl6| , |l7|| and 
in ordinary hydrodynamics |lj,|l^]. Also, the concept of 
self-affinity is not valid because of the length scale in- 
duced by ground state breakup. 

Finding the ground state of the RFIM maps exactly 
into the minimum cut - maximum flow problem of net- 
work or combinatorial optimization pOf. The use of 
such algorithms, pioneered by Ogielski |21| , has recently 
started become common as one can do exact disorder 
averages for systems governed by zero-temperature and 
energy landscape effects 22 . A related problem solvable 



with the method is the DAFF (diluted antifcrromagnet 
in a field) providing an experimental realization. 

The application of combinatorial optimization starts 
by augmenting the RFIM with two extra sites. The 
network optimization problem is defined on a graph, in 
which each edge corresponds to a site in the augmented 
RFIM. Each of the original sites is connected with one of 
the two extras, depending on the sign of the local field 
hi. The capacities of the vertices in the graph are equal 
to either 2J or 2\hi\ for couplings to the extra sites. This 
is a network flow problem, since the connections equal 
local flow constraints or capacities. The maximum flow 
between the extra sites gives the ground state energy, and 
the division to two spin states among the Ising spins is 
the minimum cut that results in the maximum flow. This 
method is exact and does not suffer from metastability 
like normal Monte Carlo or optimization with simulated 
annealing. We use an efficient push-relabel preflow-type 
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H code. The CPU-time scales as t C pu ~ N 1 - 2 , with 
N ~ L 2 increasing thus almost linearly in TV. Systems 
can be studied up to L = 1000 (N = 10 6 ). 

Figure [j] and the inset show examples of ground states 
with weak and strong disorder and with domain wall- 
enforcing or periodic boundary conditions, respectively 
PH . The random fields hi obey either a bimodal distri- 
bution {P(hi) = \[8(hi - A) + 6(hi + A)]) or a Gaussian 
one, (P(hi) = l/(\/27rA) exp[i(/ii/A) 2 ], A measures the 
standard deviation, J = 1). First we characterize the 
transition from the case of Fig. [I], a ferromagnetic ground 
state here with an imposed domain wall, to that shown 
in the inset, with a negligible magnetization. Then the 
properties of single domain walls are studied. 

We make the assumption of one single length scale, 
proportional to that at which the order vanishes (e.g., the 
magnetization becomes zero) . We measure the probabil- 
ity of a purely ferromagnetic GS, Pfm(L) = P(L, \m\ = 
1), as a function of L with fixed A. This probability 
maps for both types of disorder to the magnetization 
[m = iti(Pfm)]- The break-up length scale is defined 
with Pfm(L) = 0.5. The advantage is that the breakup 
of the ground state is visible at much smaller L than 
with other order parameters, making it possible to study 
breakup to a domain structure with L — > oo. Other 
choices could be the cluster size distribution, the spin- 
spin correlation length, magnetization, and so on. For 
example, the correlation length shows finite size effects, 
which might be partly explained below. 

Lf, for this definition is depicted for varying A in Fig. |^. 
The prediction that the 2D RFIM ground state should 
have no long-range order is based on the fact that at large 
enough scales entropy, the many possible configurations 
available should make the domain wall energy vanish || . 
Our results yield, in agreement, an exponential length 
scale 

L 6 ~cxp(A[l/A] 2 ), 

where the disorder-dependent constant A = 1.9 ±0.2 and 
2.1 ±0.2 for bimodal and Gaussian disorder, respectively. 
The definition of Lf, implies that the magnetization van- 
ishes at a larger L > Lb- The values of A are differ- 
ent from the ones obtained by finite-temperature Monte 
Carlo simulations for small L J?J. These results prove 
that the mechanism for the breakup of the GS is due to 
entropic effects. 

No ferromagnetic order exists in the ensuing domain 
structure with zero magnetization. For strong disorder 
one can show that the spin-spin correlation length is pro- 
portional to the average cluster size for both L < Lb and 
L > Lb. Here we study the disorder averaged properties 
of the largest clusters. These are found to percolate and 
thus give rise to sub-dominant (the weight of the span- 
ning cluster vanishes in the thermodynamic limit) long- 
range order. For bimodal disorder the fractal dimension 
is df = 1.90 ± 0.02 (Fig. ^ inset), very close to the ex- 
act value of standard 2D percolation 91/48. The inset 
of Fig. |2| also shows the sum over the random fields of 
the percolation clusters. This sum scales with the same 



fractal dimension 1.90. Thus the Imry-Ma argument is 
not true for the largest clusters as the global optimization 
produces domains whose magnetization is extensive. To 
summarize for weak disorder there is hidden order in the 
ground state of the RFIM structure in two dimensions. 
This is not in contradiction to the exact Aizenman-Wehr- 
theorem since m — > 0. However, it gives rise to nontriv- 
ial correlations in the structure, thus order. For stronger 
fields there is a crossover to site percolation and a non- 
percolating structure (as p c ~ 0.593 on a square lattice 
and now p — 0.5). The critical A c , below which lattice 
effects are smeared out, is (h/J) c — 2 for bimodal disor- 
der and A c = 2.1 ± 0.2 for Gaussian, respectively. The 
threshold for the Gaussian case is a rough estimate. It 
would interesting but hard to analyze this percolation 
transition in detail, since one needs L > Lb- 

Next we turn to interface scaling. Fig. || shows the 
interface width, the interface energy E, and the energy 
fluctuations AE 2 = (E 2 -E 2 ) up to L = 500-1000. E is 
obtained in the DWRG sense by subtracting the energy 
of a ground state from one with an imposed domain wall 
and identical disorder. (...) is the average over disorder. 
We take the Solid-On-Solid (SOS) -limit: in the case of 
a multiply valued interface the highest location is chosen 
from the exact interface configurations. The weight of 
overhangs is negligible for weak disorder and small sys- 
tems. In the weak-disorder regime the global roughness 
exponent is found as Q — 1.2 ± 0.05. As £ > 1 the RF 
interfaces are super-rough. This is, however, true only 
up to a length scale, below which the GS has already 
broken down (see inset, Fig. |l|). Above that scale a do- 
main wall becomes fractal, and £ = 1. There is a sharp 
transition between these two regimes, and the data for 
A = 10/9 in the inset of Fig. ||(a) has two regimes cor- 
responding to "ferromagnetic" and "disordered" ground 
states. Figure |(b) shows the DWRG result for the DW 
energy: there is a logarithmic correction to the DW en- 
ergy in the FM phase. In the paramagnetic phase the 
energy has only a remnant contribution from the bound- 
ary conditions. For the FM phase the energy fluctuation 
exponent is 6 ~ 1. The values for the exponents £ and 6 
disagree with the exponent relation 6 2 d = 2C2d — 1 |Tl] ]. 

If one studies interfaces based either on a mapping to 
the Burgers'-KPZ equation or functional RG calculations 
these depend on the small slope approximation, which is 
a problem if £ ~ 1. Fig. [| shows the statistics of interface 
fluctuations in the form of the interface step height prob- 
ability density function (PDF) /(Az^ i+ i,L). Azj,j+i is 
the height difference between two neighboring sites (zi) 
along the SOS interface (i = 1,...,L). The f(Az,L) 
show stretched exponential behavior. The PDF's are 
clearly L dependent, but only up to the breakup length 
scale for interfaces. The height differences resemble ve- 
locity gradients or energy dissipation in fluid turbulence 
and behavior in interface growth problems p^ , rf7| , p5| , p^[ 
that are governed by intermittent, rare events. The step 
height fluctuations are not restricted to the exact inter- 
faces. A SOS transfer matrix calculation (allowing for 
Az > 1) reproduces these features, and demonstrates 
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that the interpretation of Ref. |l4j is wrong since the 
true scaling behavior is super-rough at, also, low temper- 
atures. A multifractal study of the average step height 
| Ax | and the interface height-height correlation functions 
Gfc(r) = (\h(l)h(l + r)\ k ) indicates that the local inter- 
face scaling is multiafhne, e.g., Gk(r) ~ Afcr a * . For in- 
stance, ah — 0.88, . . . , 0.9 for k = 1, but for already k = 2 
as, ~ 0.66 with the exponent being a weak function of L 
at fixed A. The higher exponents a% decrease with k and 
increase with L for moderate L. Thus, there is another 
analogy between height-height correlations of the RFIM 
and velocity-velocity correlations in turbulence. Both the 
amplitude of G (Ak) and \Az\ do not self-average but 
scale with L and A. It is tempting to draw a parallel 
with Lb in here and the outer length scale of turbulence. 
In both cases the largest length scale is fixed by external 
conditions: the strength of randomness or the Reynolds 
number The correspondence is not one-to-one, how- 
ever, since the scaling properties depend on both Lf, and 
L (Fig. 4). 

In conclusion, we demonstrate the breakdown of the 
ground state, at zero temperature, in the 2D random field 
Ising model. There is however hidden, long-range order 
in the form of the spanning cluster that seems first con- 
tradictory to destroyed ferromagnetic order. This arises 
by "entropic optimization" so that the cluster magne- 
tization becomes extensive. The annihilation of order 
with increasing sample size is reflected in the properties 
of domain walls. For small systems and weak fields the 
domain walls are super-rough, with a roughness expo- 
nent that is well in excess of analytic estimates. This 
can be traced to "turbulent" rare interface fluctuations, 
but in an equilibrium system in contrast to models of ki- 
netic surface roughening or Navier-Stokes turbulence. In 
large systems the concept of an individual domain wall 
becomes ill-defined. The domain wall energy has a loga- 
rithmic correction: one should study how far this lack of 
self-affinity penetrates the GS energy landscape proper- 
ties and, perhaps, dynamical behavior. It will be inter- 
esting to see if the nonstandard interface scaling prop- 
erties persist in higher dimensions or in the presence of 
an applied field. We believe that this is so in the lat- 
ter case, though the ground states are naturally ferro- 
magnetic. This would have consequences for driven in- 
terfaces below the crossover to annealed disorder for a 
strong enough driving force |^5],|7| . 

This work has been supported by the Academy of Fin- 
land (MATRA and M. J. A.). 
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FIG. 1. Two typical examples of RFIM ground states with 
periodic and forced boundary conditions. Weak bimodal dis- 
order, A = 10/17, L = 100. Note the large jumps on the inter- 
face and the lack of overhangs. The interface is the boundary 
between the black/white domains (different spin states). In- 
set shows a strong Gaussian disorder case, A = 2. Spins that 
point "down" in the ground state are drawn white. The "up" 
spins are black or grey if in the largest (percolation) cluster. 
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FIG. 2. L;, vs. (1/A) 2 for bimodal and Gaussian disor- 
der (closed circles and open squares, respectively), calculated 
from PFM(Lt) — 0.5. A = h/J for binomial and A = 5h 
for Gaussian disorder. The inset shows the average mass of 
spanning clusters for bimodal A = 25/13 up to L — 470 (open 
diamonds) . The plot shows also the sum of the random fields 
of the sites belonging to the same clusters (closed triangles). 
The 2D percolation fractal dimension df = 91/48 is indicated 
with a line. 
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FIG. 3. (a) Scaling of the global interface width for bi- 
modal disorder. A = 2/3 (open triangles) and 3/2 (closed 
squares). The line indicates a least-squares fit with a rough- 
ness exponent £ = 1.20 ± 0.05. The inset shows the crossover 
in interface properties with increasing system size (A = 10/9). 
(b) Scaling of the energy (per length) for bimodal A = 1/3 
(open circles) and 5/12 (closed squares). The inset shows the 
scaling of energy fluctuations for A = 1/3. 




FIG. 4. Interface step probability density function for 
h/J = 1/2, L = 20, 40, 100, 200, 480. For simplicity the data 
includes only those steps that do not involve local overhangs. 



